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Animal cells can sense chemical gradients without moving, and are faced with the challenge 
of migrating towards a target despite noisy information on the target position. Here we discuss 
optimal search strategies for a chaser that moves by switching between two phases of motion (run 
and tumble), reorienting itself towards the target during tumble phases, and performing a persistent 
random walk during run phases. We show that the chaser average run time can be adjusted to 
minimize the target catching time or the spatial dispersion of the chasers. We obtain analytical 
results for the catching time and for the spatial dispersion in the limits of small and large ratios of 
run time to tumble time, and scaling laws for the optimal run times. Our findings have implications 
for optimal chemotactic strategies in animal cell migration. 


Eukaryotic cells migration does not always occur con¬ 
tinuously: bimodal motions with reorientation phases 
where cells loose their polarity and phases of polarized, 
persistent motion have been reported in unicellular or¬ 
ganisms in cell migration occurring during develop¬ 
ment in the Zebrafish and Xenopus embryo [Mg , or in 
mammalian cells |g. These two phases of motions have 
been denominated runs and tumbles [5] in analogy with 
E. Coli motion, although the corresponding migration 
mechanisms and chemical sensing of animal cells are very 
different to the ones employed by bacteria. Possibly, cel¬ 
lular run and tumble behavior reflects the necessity for 
the cell to retract motile protrusions before forming new 
ones |5]: in that case, the run time is associated to the 
protrusion lifetime (cf. Fig[^). 

Primordial germ (PG) cells in the Zebrafish embryo 
in particular have been reported to switch between two 
behavioural modes denoted “run” and “tumbles” [2]. In 
addition, PG cells move collectively in response to the 
chemokine SDF-la towards the gonad during the first 24 
hours of development [2] ■ In contrast to bacteria, after a 
tumble phase, PG cells choose a direction biased towards 
their target : animal cells are indeed capable of migra¬ 
tion towards a source through sensing of chemical gradi¬ 
ents ElIH], allowing them to bias their motion towards a 
chemical source liiin] without the need to sample the 
chemoattractant at different spatial positions. 

PG cells therefore display intermittent, directed migra¬ 
tion towards a target. We investigate here the chemotac¬ 
tic efficiency of this class of motion with a minimal model. 
We ask in particular (i) how does the time necessary to 
find the target depend on the properties of the cell mo¬ 
tion, and (ii), how does the spatial dispersion of a group 
of cells evolve in time. Both questions are relevant for 
biological processes, where cells have to move to specific 
locations and maintain their cohesion m- 

In random search processes, a chaser without infor¬ 
mation on the target location performs a random walk 
until hitting the target by chance naiig. In this con¬ 
text, the persistence length of a particle performing a 


random search can be optimised to minimise the search 
time El- In this work, we consider instead the optimal 
moving strategy for another class of motion, directed ran¬ 
dom search, where the chaser has noisy information on 
the location of the target, and needs to slow down to 
reorient. More specifically, we consider an intermittent 
moving chaser which stops for a finite time to reorient 
towards its target (e.g. due to chemotaxis), with some 
detection error. This simple model accounts for the gen¬ 
eral situation where the reorientation of the chasing agent 
takes time, either due to physical constraints or to a finite 
detection time. 

The main question is: how often should such a particle 
stop and reorient to efficiently move towards the target? 
Intuitively, short runs have the advantage of frequent re¬ 
orientation towards the target, but at the cost of frequent 
stopping, while long runs may significantly deviate the 
particle from the location of the target, due to initial ori¬ 
entation errors and possible motion of the target. We 
show here that this results in an optimal run time which 
depends on the size and distance of the target and on 
the orientation errors that the chaser makes. Similarly, 
the effective diffusion of a collection of chasers can also 
be optimized, which might be essential for collective cell 
migration towards a common target. 

We consider point-like agents (chasers) moving in two 
dimensions towards a disc-shaped target with radius S 
located at position x-r (Fig. [^. The chaser position is 
denoted X (7 and its direction of motion is given by the 
heading angle (p. The chasing agent switches between a 
run and a tumbling state. The tumbling phase lasts for 
a time tt during which the chaser does not move, after 
which the agent picks a direction uc towards the target, 
with an additional angular “detection” error rj (Fig. [^. 
The initial polar angle of the chaser at the start of the 
run phase then reads 

ipi = arg(xT - xc) + r]- (1) 

For simplicity, rj is drawn from a normal distribution with 
standard deviation e. In general, the detection error e 
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FIG. 1. (A) Run and tumble migration of animal cells: Runs 
correspond to directed migration, of polarized cells with actin- 
rich protrusions (lamellipodia) at the leading edge (green). 
Tumbles correspond to reorientation events associated with 
repolarization events. (B-D) Model schematic: A chasing par¬ 
ticle aims at a target by switching between persistent motion 
towards the target (runs) and reorientation events (tumbles). 
After reorientation, the chaser orients itself towards the target 
with an error 77 . (E) Simulation snapshot of 20 chaser trajec¬ 
tories for Tr = 8 and t < = 744.6; the red circles indicates 

chaser positions at maximal dispersion (tmax = 388.0). (F) 
Dispersion (t) versus time for different run times. The ver¬ 
tical lines indicate the average catching time Tc. Parameters 
(d,e): do = TT, do = 500, e = 0.3, S = 0.1, N = 40000. 
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FIG. 2. Gatching time Tc (a,c,e) and maximum dispersion 
<^max (b,d,f) versus r^/Tt for different orientation errors e 
(a,b), target radius S (c,d) and angular noise (e,f); in¬ 
sets show the effective diffusion coefficient DeS = f^max/tmax- 
Default parameters: N = 10^, 60 — n, e = 0.2, do = 1000, 
Tt = 1, nc = 1, S = 1, = 0. 


might depend on the distance to the target d, for instance 
due to the spatial variation of the chemoattractant. In 
this letter, we focus mainly on the simplest case of a 
constant e, and relax this assumption at the end of the 
letter. During the run phase lasting for a time tj., the 
particle is not able to further bias its motion towards the 
target, and moves with a constant speed vc and heading 
angle dynamics = ric{t), with r]c{t) a Gaussian white 
noise with variance D^. For = 0 the chaser moves 
on a straight line in the direction initially chosen ip = cpi, 
whereas for D^p > 0 it performs a persistent random walk. 

We assume here random switching between the run 
and reorientation phases to be given by a dichotomous 
Markov process. The durations of the different phases, 
tr and tt, are exponentially distributed with an average 
run time and an average stopping time Tj. 

To gain insight in the dynamics of the chasing motion, 
we performed simulations of the motion of a collection 
of chasers (Fig. [^. At time t = 0 the chaser particles 
are placed at an angle Oq and distance do from the target 
placed at the origin, xt = 0. Each chaser starts with a 
run with initial heading given by Eq. Timescales are 
set by the tumbling time Tt, velocities are normalized by 
the chaser velocity vc and distances by vcTt- 

To characterize the motion of chasers, we numerically 
evaluate the average catching time Tc and the spatial dis¬ 


persion in chaser positions, = {{xc{t))^) —{xc{t))^■ 
The dispersion exhibits different regimes over time (Fig. 
1, |15j): an initial diffusive spread, where increases 
roughly linearly, is followed by a transition to a station¬ 
ary state, where assumes a constant value. 

We then evaluated the average catching time Tc, the 
maximum dispersion during the catching process cr^axj 
and the initial effective diffusion Dgs = cr^ax/^max with 
tmax being the time of maximum dispersion. When vary¬ 
ing the ratio of run to tumble time we observe that 

Tc is minimal for an optimal value of TcjTt (Fig. [^. cr^ax 
also exhibits a minimum for zero or a finite value of the 
run time r^.. This behavior can be understood intuitively 
as follows: for short run times t^, the chaser reorients fre¬ 
quently, allowing it to follow accurately the target; the 
frequent reorientations however slow the chaser down. 
For long run times on the other end, the chaser has a 
high probability to miss the target due to orientation 
errors, leading to large displacements away from the tar¬ 
get. When varying the rotational diffusion of the chaser 
D,p, we find that for the chaser “forgets” its 

initial direction during each run and the catching time 
sharply increases, while for short run times the catch¬ 
ing time and chaser dispersion are unchanged (Fig. |^. 
Because rotational diffusion then essentially introduces a 
maximum value of above which the catching process 
becomes undirected, we focus on the limit D,p —>■ 0 in the 
following. 
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FIG. 3. Limit cases of the chaser motion, (a) For small 
enough runs and away from the target, the chaser undergoes 
a nearly one-dimensional motion, (b) For large runs and close 
to the target, the chaser moves around the target and has a 
fixed probability of hitting the target at each run. 


To understand the origin of the optimal behavior, we 
note that two regimes emerge from the analysis of nu¬ 
merical simulations depending on the ratio between the 
average run length and the distance to the target (Fig. 
[^. For VcTr |xc| (large distances and short runs), 
the chasing motion is dominated by a slow effective drift 
towards the target. For VcTr > Ixd (small distances and 
large runs), the chasers approach the vicinity of the tar¬ 
get and the catching process is primarily dominated by 
the probability of hitting the target within a single run. 
The total catching time is then the sum of the time to 
approach the target (first regime), and the time to catch 
it (second regime). A visualization of the catching pro¬ 
cess for different average run times is provided in the 
Supplementary Material Ha- 

Based on these observations, we start by considering 
the limit of short runs with chasers positioned far from 
the target with the initial distance do 3> S. We assume 
here, without loss of generality due to rotational symme¬ 
try, that the chaser initial position is along the x axis, 
Xc(f = 0) = —doe^- The chaser then follows a nearly 
one-dimensional motion along the x-direction, such that 
the desired chaser heading angle is small ((^ « 0). The 
system can be described in terms of the spatial probabil¬ 
ity density functions (PDFs) for chaser particles in the 
run and tumble phase Pr(x, (p) and pt(x) [TB], where the 
probability distribution is conditioned to the initial posi¬ 
tion and angle (xq, po)- The general evolution equations 
for the PDFs in the frame of the target read 


dtPr{yi,p) 


dtPti^) 


Ur.Vprix, (p) - —Pr{x, p) + —T(p)pt{x) 
Tr Tt 


-Pt(x) -b — / 
'^r Jo 


27T 


dpPr{x,p) 


( 2 ) 

(3) 


where = vc'^cip) is the drift velocity in the run 
phase, T{p) = ^^/j^exp is the probability of re¬ 

orientation in the p direction after a tumble, V denotes 


the gradient operator acting on the x dependency of 
Pr{x, (/?), and the time dependency of Pr and pt is implic¬ 
itly implied. By deriving moment equations from Eqs. 
2|3 (Supp. Mat), we find that the average chaser veloc¬ 
ity towards the target relaxes in a time TVTt/(rr -b Tt) to 
the stationary value 


(v) = 


VcTrC 


Tr + Tt 


(4) 


while the chaser dispersion fT^(t) = Degt is predicted 
to increase linearly in time with the effective diffusion 
coefficient 


DeS 


‘^Vc'^r 

Tr + Tt 


1-(1 


(Tr + Tt)"^ 



(5) 


For e = 0, in timt limit, the diffusion is 

only due to the stochasticity in run and tumble times, and 
for large run time the diffusion constant decreases with 
Tr, since each additional run introduces dispersion by 
generating different run lengths between different chasers 
(a zero run time also minimizes the diffusion, but the 
catching time then diverges). For e 7 ^ 0, there is an opti¬ 
mal run time minimizing the diffusion coefficient, given 
for e —)■ 0 by Tr/rt = 1 /e: large runs indeed increase dis¬ 
persion, because they tend to amplify initial errors in the 
heading angle. 

The average approach time to the target is given in 
that drift-dominated limit by 


T(i) = 


do-S 


VC 


Tr+Tt 


( 6 ) 


which agrees with numerical simulations (Fig. SI). 
Therefore, in the limit of small runs, increasing the run 
time is always favorable. 

We now discuss the opposite limit when the target is 
sufficiently close to the chaser to be reached within a sin¬ 
gle run (Fig. [^). We assume that the probability distri¬ 
bution of the chaser position at tumble, pt, has relaxed 
to a steady state. The steady-state distribution is then 
isotropic. At lowest order in the orientation error e, the 
chaser-target distance R = |xc| is distributed according 
to an exponential distribution, pt{R) = -d—e~ , with 
the corresponding standard deviation of chaser position 
CTstat = TrVc- To evaluate the catching time, we note that 
the probability of hitting the target in one run, after each 
tumble, is given by 

pOO 

Phit — I dRpt{R)pdiT{R)PTeach{R) ( 7 ) 
Jr=o 

where pdir (R) is the probability of choosing a direction 
towards the target, and PreachiR) is the probability of 
performing a sufficiently long run to actually reach the 
target. We start by evaluatingpdir- The chaser is heading 
towards the target provided that the detection error is 
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FIG. 4. Comparison of simulation results and theory: (a) 
Tc and analytical prediction Tc = (b) (Tmax with 

DesTc for the “quasi-lD” approach and o-gtat = (TrVc)^ for 
the “single run” encounter (right). Parameters: N = 10^, 
e — 0.1, S' = 1, do = 500. Scaling of the optimal run time t* 
minimizing Tc, (c) with the angular error e, (d) the average 
tumble time rt, (e) the target size S and (f) the initial distance 
d. The solid lines indicate the predicted dependence from 
theory. Simulation parameters: N = 10"^, Vt = 0.1, and 
unless parameters are varied, rt = 1, vc = 1, e = 0.2, S = 1, 
do = 500 and 9o = tt. 


sufficiently small, \r]\ < As a result, the probability of 
choosing the right direction by a favorable choice of error 
r] is given by 


Pdir{R) = 2 




( 8 ) 


To obtain Preach(.R), we note that (1) at lowest order in 
the run time has to to be larger than fcatch = R/vc 
to reach the target, and (2) run times are taken out of an 
exponential distribution; therefore, Preach = . 

Overall, we then obtain phit — 

M{x) = dye~y'&Yi{x/y) [TT], and the average catch¬ 
ing time from the vicinity of the target reads 


^(2) 


2(Tr +Tt) 
'■ TrVc^ ' 


(9) 


Therefore, in the limit of large runs, decreasing the run 
time is favorable, an opposite behavior to the short run 
time limit discussed before. 

The total time to reach the target, Tc, is the sum of 
the time required to approach the target (Eq. and the 
time to catch it in its vicinity (Eq. |^, Tc = Th ^ -|- Tc^\ 
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FIG. 5. Catching time Tc and maximum dispersion 
versus ei. Simulation parameters: N = 10"*, Tt = 1, uc = 1, 
£0 = 0.1, S' = 1, do = 300 and 9o = tt. 


The resulting analytic prediction accounts very well for 
the catching time obtained by simulations (Fig. |^). To 
estimate the optimal run time, we look for the minimum 
of the total catching time Tc. We find that for interme¬ 
diate orientation error e 1, the optimal run 

^dVaTt 

time T* scales as m 



( 10 ) 


The condition can be rewritten as an optimal run length 
T*vc ~ {dSrtVc/d)^, scaling with the geometric mean of 
the three lengths of the problem. The predicted scaling 
is verified by numerical simulations (Fig. |^-f). The 
intermittent directed chaser has an optimal run length 
strikingly different from the optimal persistence length of 
a random searcher, which was found to be of the order of 
the system size [14]. It further implies that performing 
short length runs is better for small target sizes, when 
the distance to the target is close, and for short tumbles. 
The optimal run length is also crucially dependent on 
the orientation error: quite naturally, smaller orientation 
errors allow to perform larger runs. 

In the minimal model discussed so far, we have as¬ 
sumed that the reorientation error e is independent of 
the distance to the target d{t) = |xc(t)|. For directional 
sensing however, the detection accuracy of the target di¬ 
rection is likely to decrease with d{t). To account for this, 
we numerically investigated an extension of our model 
by assuming that an arbitrary error function e{d) can 
be approximated by a linear function e{d) « eg + eid 
with €i > 0. Simulations results for a finite ei show the 
same qualitative behavior as the minimal model (Fig. 
[^; an optimal run time minimizing the catching time 
still exists. When increasing ei, the optimal run time 
minimizing the dispersion moves towards zero, similar to 
increasing the overall constant error e (Fig. 

We have also considered the situation where the target 
moves balistically with a finite, but small, speed vt < vc, 
and found that the same qualitative findings hold m- 
A crucial difference for a finite target speed vt is the 
existence of a regime where the catching time Tc diverges. 
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occurring above a threshold detection error e or for short 
run time Tr [15] . 


The robustness of our qualitative results to these mod- 
ihcations suggests that the minimal model introduced 
here exhibits physical properties which will be preserved 
for more realistic description of cells moving with inter¬ 
mittent motion. Developmental processes require cell mi¬ 
gration occurring at the scale of the organism HHj: it 
would be interesting to experimentally test whether mi¬ 
grating cells operate near the optimum value predicted 
by Eq. 10 in this context. Recently, Minina et al. [3] 
have analyzed in detail guided migration of progenitor 
cell in vivo. Based on their experimental results they sug¬ 
gest that progenitor cells might adaptively decrease their 
run length by downregulating receptor signaling close to 
the source of the relevant chemokine. A disruption of 
this control mechanism for run length was shown to hin¬ 
der precise arrival of cells at the target. These empirical 
findings are confirmed by our theoretical work, in partic¬ 
ular by the predicted decrease in optimal run length with 
decreasing distance to the target (Fig. [^). This indeed 
suggests that PG cells have evolved complex chemotactic 
strategies for precise spatial and temporal arrival at tar¬ 
get sites, to ensure successful development of the embryo. 
We expect that further research combining new quantita¬ 
tive experiments and more elaborate theoretical models 
will provide important insights into the mechanism and 
function of directed cell migration in living systems. 


From a more general point of view with respect to the¬ 
ory of search processes, our simple but generic model can 
be viewed as an investigation of optimal directed search, 
which is fundamentally different from the random search 
processes predominantly studied up to date. 
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